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SUMMARY 


A further examination of previously reported unsteady aerodynamic data 
on a tunnel spanning wing (both swept and unswept), obtained in the United 
Technologies Research Center (UTRC) Main Wind Tunnel, was performed. In the 
original study (NASA CR-3092) the analysis of these data was focused on the 
behavior of the integrated force components of the NACA 0012 airfoil section 
tested (lift, drag, and pitching moment). The objective of the present 
investigation was to study the unsteady chordwise force response on the air- 
foil surface and to examine its sensitivity to the various system para- 
meters . 

The main body of this data analysis was carried out by analyzing the 
propagation speed of pressure disturbances along the chord, and by studying 
the behavior of the unsteady part of the chordwise pressure distribution at 
various points of the airfoil pitching cycle. It was found that Mach number 
effects dominate the approach to stall and the inception of stall, both 
static and dynamic; the stall angle decreases as the Mach number increases. 

It was also found that sweep dominates the load behavior within the stall 
regime; large phase differences between unswept and swept responses, that do 
not exist at low lift coefficient, appear once the stall boundary has been 
penetrated . 

Conventional pressure normalization procedures were reviewed and it was 
found that the cosine law for sweep angle is not applicable at angles of 
attack above the stall angle, both for steady and unsteady conditions. It 
was also found that the reduced frequency parameter traditionally used in 
unsteady aerodynamic analysis is not a reliable indicator of the unsteady 
aerodynamic response in the high angle of attack regime. 

An additional analysis was performed to resolve a troublesome aspect of 
the previously reported results. It was found that the relative displacement 
of the lift and moment hysteresis loop response (excluding nonoscillatory 
terms) during upstroke between unswept and swept wings is a real effect, and 
is almost entirely described by the phase and amplitude differences in the 
first harmonic components. 
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INTRODUCTION ■ 


The need for swept wing unsteady aerodynamic data into the dynamic stall 
regime was addressed in Ref. 1 wherein test results were reported from exper- 
iments performed using a NACA 0012 airfoil section tunnel-spanning wing (TSW) 
oscillating in pitch in both the unswept and 30 deg swept back configura- 
tions. The main motivation for obtaining such data is for use by the 
helicopter designer who must account for a widely varying aerodynamic 
environment on the rotor blade as it traverses the azimuth plane. In addi- 
tion to this immediate engineering application, these results provide the 
opportunity to examine the effects of sweep at moderate pitching amplitudes, 
and thus permit the evaluation of techniques currently in use for estimating 
the unsteady swept wing aerodynamic forces from the existing data banks for 
zero sweep angle. 

The results that were reported in Ref. 1 focused on (1) the unsteady 
integrated load components (lift, drag, and moment) computed from unsteady 
chordwise pressure distributions for both stalled and unstalled conditions, 
and (2) the surface pressure fluctuations obtained in the stall regime from 
companion steady airfoil experiments. Several findings were reported, of 
which the most notable is that sweep tends to delay the onset of dynamic 
stall. A summary of the main results can be found on page 40 therein. A 
companion data report to Ref. 1 was also published (Ref. 2) containing a 
complete set of unsteady normal force, chord force, lift force, drag force, 
and moment responses in the form of hysteresis loops. The data report was 
arranged to emphasize the effect of sweep on the results. A review of 
related work in this field can be found in the introduction to Ref. 1 and in 
the associated bibliography. 

The present report utilizes the same data base that was used to produce 
the results in Refs. 1 and 2 and is essentially a continuation of that data 
analysis. Here, the emphasis is on the unsteady chordwise pressure data and 
their sensitivity to the various test parameters: wing sweep angle. A; 

pitching amplitude, a; mean angle of attack, a^; pitching frequency, f; and 
Mach number, M c = cos A. Two parametric studies are presented herein: 

(1) the behavior of the chordwise propagation of pressure waves associated 
with the dynamic stall process, and (2) the behavior of the unsteady surface 
pressures, both individually as time histories over the whole motion cycle, 
and collectively as chordwise distributions at specific points of the motion 
cycle. 

An accompanying data report (Ref. 3) contains listings of the first ten 
harmonic components of the individual chordwise pressure responses, and of 
the integrated load results. The data tables in Ref. 3 can be used to repro- 
duce the oscillatory components of the individual pressure time histories, 


2 



and the complete normal force, chord force, lift force, drag force, and 
moment coefficient time histories. (Here, the term "complete" means 
nonoscillatory plus oscillatory parts, and it is understood that the chord 
and drag forces were computed by integration of pressure distributions and 
therefore do not include the effects of viscous drag.) Hysteresis plots of 
the ingegrated force coefficient results (versus angle of attack) were 
reported in Ref. 2. 

The present document begins with a review of the data base and the 
approach to the normalization of the data. This is followed by a description 
of the data retrieval process that was used to expedite the performance of 
the parametric analysis. A discussion of the chordwise wave propagation 
phenomenon and its sensitivity to the test parameters is then presented. A 
review of the effectiveness of the data normalization procedure is presented 
next. The report then proceeds with the analysis of the chordwise response 
behavior to the various test parameters, including a brief examination of the 
individual pressure time histories very near the leading edge. Finally, a 
troublesome aspect of the hysteresis loop comparisons of Ref. 1 is revisited 
and resolved. 
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LIST OF SYMBOLS 




■'D 


"M 


'N 


c p ( x,t ) 


f 

k c 

M c 

M» 

m 

p(x,t) 

q 

T 


dimensionless pitch rate, Eq. (14) 

dimensionless pitch rate at vortex inception angle 

sound speed 

airfoil chord 

chord force coefficient 

drag force coefficient, Eq. (5) 

lift force coefficient, Eq. (4) 

moment coefficient 

normal force coefficient 

pressure coefficient distribution, Eq. (3) 
airfoil pitching frequency, Hz 
chordwise reduced frequency, Eq. (2) 
chordwise Mach number, Eq. (1) 
free stream Mach number 
slope, Eq. (11) 

o 

chordwise pressure distribution, newton/m 
dynamic pressure based on chordwise velocity, newton/m^ 
period, sec 
time, sec 

vortex inception time, sec 
chordwise velocity, ra/sec 
free stream velocity, m/sec 
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LIST OF SYMBOLS (Cont'd) 





a 

a 

a 

a 

a 

a 

e 


M 

s 

ss 

V 

VO 


A 


T 


T 


O 


X 

0) 


average dimensionless wave speed, Eq. (8) 

average dimensionless wave speed between stations i and j 
airfoil angle of attack, deg, Eq. (6) 
mean angle of attack, deg 

dynamic stall inception angle, deg, Eq. (19) 

steady state stall angle, deg, Eq. (19) 

vortex inception angle, Eq. (10) 

constant used in Eq. (11) 

data scatter parameter, Eq. (15) 

sweep angle, deg 

dimensionless time based on oscillation period, also 
fraction of period 

dimensionless vortex inception time based on oscillation period 

phase lead, deg 

distance from airfoil nose 

circular frequency, rad/sec 
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DESCRIPTION OF THE DATA BASE 


The data base for the present study was obtained from aerodynamic exper- 
iments performed on an oscillating tunnel-spanning wing (TSW) in both unswept 
and 30 deg swept configurations in the United Technologies Research Center 
(UTRC) Main Wind Tunnel. The chordwise profile taken normal to the wing 
leading edge is a NACA 0012 airfoil. A detailed description of the test 
program can be found in Ref. 1. Although the data base and normalization 
procedure are also fully described therein, the portions that are pertinent 
to the present study are repeated here for completeness. Of primary interest 
in this investigation is the unsteady chordwise pressure distribution. The 
unsteady data acquisition points that were used along the chord are schem- 
atically depicted in Fig. 1. 

The experiments of Ref. 1 were performed for the matrix of test para- 
meter settings shown in Table 1. These parameters were the wing sweep angle, 
A, the amplitude of pitching motion, a, the mean angle of attack, a^, the 
pitching frequency, f, and the approach Mach number, 1^. The effect of 
sweep was accounted for in the normalization process by referring all of the 
data to the component of Mach number, M c , normal to the wing span. This 
quantity is related to the freestreara value, ^ , by the formula 

M c = ^ cos A (1) 

where subscript c denotes the chordwise direction. Data were obtained for 
two values of M c (= 0.30 and 0.40) for each value of A (=0 deg and 30 
deg). This approach required testing the 30 deg swept configuration at a 
higher freestream velocity as shown in Table 1. Unsteady pressure data were 
obtained for two values of a (8 deg and 10 deg) at four mean angle of attack 
settings (0, 9, 12 and 15 deg) about the quarterchord axis. The angles were 
measured in a plane perpendicular to the span axis. The choice of pitching 
frequencies used during the test program was predetermined according to Table 
2 which lists the nominal values of reduced frequency defined by 


k 


c 


_ TT cf 
^»^C 


( 2 ) 


where is the speed of sound. 
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The effects of sweep on the pressure data were accounted for in Ref. 1 
by normalizing the data with respect to the velocity component normal to the 
span to obtain the following expression for dimensionless pressure 


C p (x»t) = - p(x,t)/qcos 2 A (3) 

where x is the dimensionless position along the chord and t is time. This 
nondimensionalizat ion scheme is based on a steady-state sweep analysis which 
assumes the flow to be potential and the pressure distribution on the wing to 
be determined entirely by the magnitude of the normal velocity component, 
cos A. Note also that in this report, the pressure coefficient is 
defined as the negative of the normalized pressure. This is convenient for 
relating positive integrated lift to a positive value of Cp. 
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DATA PROCESSING AND ANALYSIS 


Introductory Comments 

The analysis of unsteady chordwise pressures and their sensitivity to 
the various system parameters by conventional methods is a formidable under- 
taking for a data base the size of the one shown in Table 1. This is true 
mainly because an important aspect in the study of chordwise pressure 
behavior involves a comparative study of the pressure distribution at various 
points of the airfoil pitching motion cycle. This approach requires substan- 
tial data manipulat ion if a reasonable number of representat ive cases and 
comparisons are to be analyzed in any significant detail. (It is noted that 
the chordwise pressure in the present study consists of 21 channels of data, 
including 13 upper surface and 8 lower surface locations.) 

The approach that was used to expedite the present analysis consisted of 
(1) compacting the processed data into a data file where each record (test 
case) is addressable either by record number or by the corresponding test 
parameters (A, a, a^, f, and M c ) that are unique to the record, and (2) 
tailoring a data management program where specific channels from any combina- 
tion of records can be retrieved for numerical and/or graphical analysis by 
means of a series of built-in computer queries. 


Data Compaction and Storage 

The direct access data file that was used in the present analysis was 
generated in two steps: (1) compacting the unswept and swept data bases into 

a common data file, and (2) creating a file directory for interfacing the 
data file with the data retrieval software (described later in this 
section) . 

Data compaction was achieved for the unswept and swept data bases by 
cycle-averaging the individual pressure responses of each case over 5 cycles 
of the airfoil pitching motion and then reducing the resulting time histories 
to 24 harmonic components (including the zeroth harmonic term). The 
fundamental frequency used in the Fourier analysis process was the pitching 
frequency of the airfoil. The unswept and swept wing data bases were then 
combined into a common compacted data file. 

The normal force, chord force, and moment responses were then directly 
computed from separate chordwise integrations (by trapezoidal rule) of the 
in-phase and out-of-phase portions (relative to the airfoil motion) of each 
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harmonic component. The lift and drag forces were then computed from the 
normal force and chord force using the following formulas: 


cos a + C £ sin a (4) 

C D = C N sin a - C c cos a (5) 

where a is given by 

a - otjj + a sin cat. (6) 

A schematic depicting the force vectors is shown in Fig. 2. By inspection of 
these equations, it is seen that the harmonic coefficients of and Cp 

cannot be computed directly from the harmonic coefficients of and C c 

because of the nonlinear appearance of a within the trigonometric terms. 
Therefore, an intermediate computatin of the Cp and Cp responses was made 
by substituting time history reproductions of C^, C c , and a into Eqs. (4) 
through (6). The resulting Cp and Cp time histories were then Fourier 
analyzed, and together with C^, C c , and C^, were added to the compacted 

data file by extending the number of channels of each test case from 22 (a 

and 21 chordwise pressures) to 27. 

The final step in completing the contents of the compacted file was 
adding to each case the steady state force and moment coefficients 
corresponding to Cjj, C c , C^, Cp, and Cp as well as the actual test parameter 
values (A, a, a^, f, and M c ) measured during the test. 

The interface of the data file with data retrieval software was achieved 
by installing a file directory header into the data file. The file directory 
was indexed to match the actual record sequence of the file and was set up 
with nominal values of A, a, a^, f, and M c corresponding to the actual test 
parameter values of each record. The objective here is that all requests to 
the data file are processed by means of the file directory. Data retrieval 
is thus achieved without having to directly search through the individual 
records of the file. This is an especially desirable approach when records 
are being retrieved out of sequence from the file. 


Data Analysis Techniques 

The objective of the data analysis in this report is to examine and 
evaluate the effect of the test parameters (A, a, a^, f, and M c ) on the 
unsteady pressure response over the surface of the blade. To accomplish 
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this, the authors found it necessary to use several graphical representations 
while concentrating on the various characteristic features of the data. It 
is appropriate, therefore, that the descriptions of each figure type be 
consolidated in the present section so that subsequent sections can be 
devoted primarily to the discussion of the physical phenomena being 
depicted . 

Much of this report deals with changes in unsteady blade pressure 
associated with changes in certain of the test parameters. Recall that in 
this report, pressure coefficient is defined as the negative of the 
normalized pressure at any station, 

C p <X,t) - -p(x,T )/qcos 2 A . (7) 


In some figures the pressure coefficient is plotted versus dimensionless 
chord position, x> while in others the independent variable is dimensionless 
time, T. For example, Fig. 3 represents the pressure coefficient plotted 
versus x at the instant that a = (upper panel), and at the peak 
incidence angle, a = a M + a (lower panel). The direct comparison to be made 
here between the solid and dashed lines is between two chordwise Mach 
numbers. (Other plots will show the effect of other parameters in the same 
manner.) A similar chordwise pressure plot, shown in Fig. 4, depicts the 
first harmonic pressure amplitude (upper panel) and phase angle (lower panel) 
taken from the Fourier representation of the pressure time history. Here the 
direct comparison is between the two sets of symbols representing the two 
sweep angles tested. 

A third type of chordwise pressure plot is shown in Fig. 5. Here, the 
format of the individual panels is similar to that of Fig. 3, but the intent 
of the figure as a whole is significantly different. Now the solid and 
dashed lines denote upper and lower surface pressures, and the major compar- 
ison is made between the two columns. (In this case, the effect of sweep is 
depicted.) A second parametric comparison can be made in Fig. 5 by viewing 
successive pairs of panels which, in this case, are for three instantaneous 
angles of attack, a. Both the numerical value of a and the relative position 
of the blade during the motion cycle are shown between each pair. 

The last type of chordwise pressure plot shown in Fig. 6 differs from 
Fig. 5 only in the orientation of the individual plots. Here, for example, 
the upper row curves were obtained at a mean angle of attack of = 9 deg, 
and the lower row at a M = 12 deg. Each column then represents a specific 
angular displacement from the mean (i.e., 7, 7.5, and 8 deg) as the peak 
angle of attack is approached and attained during the upstroke. 

Individual pressure time histories are also used in this report. A 
typical example is shown in Fig. 7. In every case the upper row of panels 
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will represent the motion of the blade about the mean angle of attack, and 
each successive row will show the pressure time history for the specific 
chordwise measurement stations denoted by x- As in Fig. 5, the comparisons 
in Fig. 7 are made between solid and dashed lines within each panel, and 
between adjacent columns. 


In Ref. 4 it was shown that plots of unsteady pressure versus both time 
and chordwise position could be analyzed to yield quantitative information on 
the chordwise propagation of pressure waves. That technique is applied here- 
in. Plots of pressure time histories were produced at several chordwise 
positions, and in the example of Fig. 8, the data are restricted to the first 
15 percent of the chord for clarity. In this figure the vertical axis 
represents both the unsteady pressure coefficient level relative to steady 
state Cleft scale) and the chordwise location of the pressure measurement 
(right scale). The result is a pseudo three-dimensional plot of the 
Cp(x» T ) surface, referred to as a pressure carpet plot. To interpret these 
figures, refer each curve to its zero level at the appropriate X”tic mark on 
the right scale, and use the left scale as a gage to measure the pressure 
relative to its zero level. The precipitous drop in Cp in the vicinity of 
t= 0.17 to 0.2 is customarily associated with dynamic stall and is usually 
ascribed to the formation and chordwise passage of a stall cell or vortex. 


For this carpet plot a set of constant pressure contours (isobars) was 
constructed in the x i T plane as shown in Fig. 9. The two columns of numbers 
on the right represent the percent of the pressure range from minimum to 
maximum and the associated value of unsteady pressure coefficient relative to 
steady state (for example, a value of Cp - -3.0 is 4 percent of the range 
above minimum Cp, and a value of Cp = 6.0 is 74 percent of the range 
above minimum). The numbers that are appended to each contour are the 
corresponding values of Cp. As shown in Ref. 4, the direction of the 
family of contour extremes identifies the direction of wave propagation, and 
a simple estimate of the slope yields the wave speed. 

The sensitivity of the contours to relatively small timewise discrep- 
ancies in data from adjacent measuring stations revealed a minor error in the 
data set, which is described in the Appendix. The error manifested itself in 
small time or phase displacements of adjacent signals, and is occasionally 
noticeable as a zigzagging ripple in the chordwise pressure distribution, as, 
for example, in the lower right panel of Fig. 5, or in several panels of 
Fig. 6. 

As noted in the Appendix the entire data set was uniformly corrected for 
the computations made in the wave speed analysis, but the corrections were 
confined to that section only. The remainder of the analysis, including all 
of the tabulated data in Ref. 3, remained uncorrected, to preserve its 
commonality with previous results from Refs. 1 and 2. 
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Normalization Procedure 


In the present analysis normalization was achieved by dividing the 
negative of the measured unsteady pressures by the dynamic pressure, 

0.5 pV^, and by multiplying the circular frequency of airfoil pitching 
oscillation by c/2V c . Here, V c is the chordwise component of the 
freestream velocity ( V cos A) and c is the airfoil chord. 
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PRESSURE WAVE STUDY 


Wave Speed Analysis 

The use of pressure coefficient carpet and contour plots for wave speed 
analysis was introduced in the previous section, using Figs. 8 and 9 for the 
forward 15 percent portion of tile airfoil as examples. The corresponding 
carpet and contour plots for the entire chord are presented in Figs. 10 and 
11. (In this particular case, A = 0 deg, a = 8 deg, a M = 9 deg, f = 10 Hz, 

M c = 0.3 and k c = 0.125.) In both contour plots only positive values of 
Cp (cf. Eq. 7) have been included, to focus on the region of maximum normal 
force and hence on the stall wave propagation. 

Figure 12 is a contour map of the pressure surface in Fig. 8, focused on 
the first four tenths of the period, 0 < x < 0.4. In this figure, the 
connection between the contour extremum at the leading edge (xq = 0.004) with 
the extremum at the 15 percent chord location (x 2 = 0.149) by the dashed line 
labeled "overall 15%" is used for estimating the wave speed over the forward 
15 percent of the blade chord. It is seen that this set of contour extremes 
moves toward increasing time as the chord is traversed from front to rear. 

The dimensionless wave speed was shown in Ref. 4 to be related to the ratio 
of chord extent to elapsed time by the formula 

V W = (8) 

and with t = tT = 2ttt/u), and k c = cw/2V c , this becomes 


W 



(9) 


A careful study of the contours in Fig. 12 indicates that the locus of 
contour extremes is not linear but can be represented by straight line 
segments of different slopes. A two-segment estimate for this case is 
illustrated in Fig. 13. The end points of the first segment are denoted by 
(X 0 , x Q ) and (x! , Xj ) and those of the second segment by <x 1 , x 1 ) and (x 2 , 
x 2 ) . Hence, there are actually two distinct regions within the first 15 
percent of the chord that have different wave speeds. For this particular 
example, these values are 

(X , X » X ) = (0.004, 0.028, 0.149) 

0 12 

(x , t , T ) = (0.123, 0.155, 0.187) 

0 12 
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and with k = 0.125, application of the wave speed Eq. (9) yields 




02 



0.035 (first segment) 


0.146 (second segment) 


0.090 (overall, 15%) 


where the wave speed subscripts match the subscripts of the initial and end- 
point coordinates of the segment in question. Furthermore, Fig. 14 for the 
entire chord shows that another discontinuity can be defined at the 15 
percent chord station, separating the forward region from the aft region of 
the chord, implying additional speed changes along the chord. This last 
coordinate point is given by 

(X , x ) = (0.971, 0.284) 

3 3 

and two additional wave speeds can be defined. 



fr _ O U 

V w TT T - T 

03 3 0 


For convenience, the individual 
collected below. 


V w = 0.035 
01 

v w = 0.146 
12 

V w = 0.326 
23 


: 0.326 (third segment) 

0.239 (overall, 100%) 
segment wave speeds for this example 

V w = 0.090 
02 

V w = 0.239 
03 


are 


It is seen that the wave is slowest in the region of high pressure gradient 
near the leading edge (Vy ). This is the region characterized by the dense 
array of constant level contour lines in Fig. 13. Conversely, the wave 
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travels more rapidly over the second segment (Vy ) where the chordwise 
pressure gradient is lower, and attains its highk&t speed over the last 
segment ) where pressure gradient is minimal. However, without 

additional ^idence, it is impossible to establish a cause-and-ef f ect 
relationship between the pressure gradient and the speed variation. 
Alternatively, starting inertia in the vortex rollup associated with dynamic 
stall could be the reason for wave acceleration. Whatever the cause, a 
consistent observation has been made (as will be shown presently), that the 
wave speed seems to be inversely related to the unstalled chordwise pressure 
gradient . 


Before proceeding with a detailed examination of parametric trends, 
several values of wave speed for the unswept wing tests are compared in Fig. 
15 with those obtained in Ref. 4. Overall averaged wave speed for the entire 
chord (Vrj ), and local averaged wave speed for 15 percent of the chord 


03 

(Vr. ), are plotted versus reduced frequency in the upper and lower panels of 
Wqo # # 

Fig. 15, respectively. The original data are denoted by the open symbols, 
and the current data for a = 8 deg and M c = 0.3 by the solid symbols. 

It is seen that the current results are consistent with previous calcula- 
tions . 


A study was performed to examine the effect of each of the test para- 
meters on wave speed and the results are shown in Figs. 16 through 18. In 
each figure the overall wave speed, V.. , is plotted versus reduced fre- 

quency, k c . Mean angle of attack variation is always depicted by the three 
symbols used: circle, triangle, and square. In each figure there are four 

panels, with one parameter held constant from left panel to right panel, a 
second parameter held constant from top panel to bottom panel, and the varia- 
tion of the parameter of primary interest for each figure highlighted by the 
contrast between open and solid symbols. In Fig. 16 the effect of amplitude 
variation is examined by noting that the open symbols represent a - 8 deg 
and the solid symbols represent a = 10 deg. This figure shows little or 
no effect of amplitude on overall wave speed. Although it was originally 
suggested in Fig. 20 of Ref. 4 that increasing amplitude appeared to increase 
the value of Vrj at any k_ , the present results do not corroborate this 

• 0 3 . c a # 

earlier indication. It is noteworthy that the data in Ref. 4 were obtained 
for a variety of Mach number values and of profile shapes, whereas the pre- 
sent results are taken from a wholly self-consistent data base. Also, this 
figure and subsequent figures show that there is little or no effect of mean 
angle of attack on wave speed, which agrees with the observations made in 
Ref. 4 . 

Mach number is the parameter of primary interest in Fig. 17, and its 
effect, although small, seems to be consistent. For A = 0 deg (top panels) 


15 


the higher Mach number (solid symbols) yields a slightly smaller value of 
than M = 0.3, although the two sets of data appear to converge as k 


. W 03 

increases. For A = 30 deg (bottom panels) there is less distinction between 
the two sets of Mach number data, and it could be argued that the effect of 
Mach number is minimal. 


By far, the parameter that has the greatest effect on Vy is the sweep 
angle, as shown in Fig. 18. In all four panels the overall wave speed for 
the unswept wing (open symbols) is consistently and significantly greater 
than that for the swept wing (solid symbols) by a factor of 1.5 to 2.0. This 
phenomenon of lower chordwise wave speed in swept flow has not been predicted 
by any previous analysis nor observed in any previous experiment. 


In an effort to examine these trends more closely, wave speeds for the 
several segments (defined earlier) were computed and are displayed in Figs. 
19 and 20 for M c = 0.3 and 0.4, respectively (at the same value a = 8 
deg). The three left hand panels in each figure are the three local wave 
speeds, V, 


W 


and V, 


proceeding aft along the chord. The two right 
«0i- "12' w 23" — 

hand panels are tne two average wave speeds over the first 15 percent 


W 


i TJ ClLLVa Vtt , Ul. J.LULU U1C ICOU1L9 1UL V rj • Virtually all of the 

w a i ^ Wi 2 0 # Wq2 

ifrect is perceived to occur in the aft portion of the chord (lower 


Wq2 

of the chord, and v w„, over the entire chord. (The latter are repeated in 
these two figures from the top panels of Fig. 18.) It is clearly shown here 
that there is substantially no effect of sweep angle on the local wave speed 
over the forward 15 percent of the chord, either by examination of the panels 
for V v , and V tt , or from the results for V T 
sweep el 

left panel of both Figs. 19 and 20), with factors of as much as 2 observed in 
swept wing wave speed versus that for the unswept wing. This result is 
representative of all cases examined. If the pressure wave is associated 
with the formation and chordwise transport of a vortex, these results would 
imply that the vortex motion along the chord is independent of sweep over the 
first 15 percent region of the chord, and is strongly dependent on sweep 
thereafter . 


Finally, an analysis of the chordwise variation of local wave speed is 
provided by Fig. 21 for the case with M c = 0.4, a = 8 deg and a M = 15 deg. 

In this figure the local wave speed is plotted versus dimensionless chord 
position, x» for several values of k. The unswept results are plotted in the 
left panels and the swept results in the right. Linear scales are used in 
the upper panels, and two effects clearly emerge. The first is the general 
increase of with k at any given value of x (although there is some 
scatter) for both values of A. The second is the continual increase with x 
Of v w for A = 0 deg, and the much smaller increase with x of aft of x = 

0.1 for A = 30 deg. In an effort to spread out the many points plotted at 
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small values of x> the same data are repeated in the semi-log plots in the 
lower two panels. While there is no unusual revelation afforded by these 
semi-log plots, they do appear to resolve the cluttered points into an 
orderly chordwise progression. 


Vortex Inception Angle Variation 

In pursuing further the interpretation of the phenomenon just described 
as the formation of a vortex and its subsequent chordwise transport, the 
inception time can be studied and related to the vortex inception angle, 
a v . A direct computation of a v , can be made by substituting the dimen- 
sionless inception time, t q = (Dt 0 /2iT (cf. Fig. 13 or 14), into Eq. (6) 
for instantaneous angular position, 

ct v - ctj^ + a sin 2 ttt 0 (10) 

The results are shown in Figs. 22 and 23 as plots of vortex inception angle 
versus reduced frequency. As before, the use of circles, triangles, and 
squares denote =9, 12, and 15 deg, and open and solid symbols highlight 
the comparison of sweep angle effects in Fig. 22 and of Mach number effects 
in Fig. 23. These plots also contain information on the steady state stall 
angle (triangular symbols pointing to the appropriate values on the ordinate) 
and on the best linear fit to the data. The latter were obtained by a least 
squares procedure for each data group, and the former by estimating the 
steady state angle of attack at which dCM/da = 0 in Figs. 14 through 17 in 
Ref. 1. 


In both figures the vortex inception angle, oty, varies linearly with 
k c , as indicated by the good correlation of the linear least squares line 


a 


v 


a 


vo 


+ mk c 


( 11 ) 


with the measured data, and the narrow scatter range for each data set. (In 
this equation ay Q is the zero frequency intercept and m is the slope.) 

Figure 22 shows that sweep has a tendency to delay vortex inception, but this 
delay is less than two degrees overall, and less than one degree in the 
majority of cases examined. In contrast to this result, Fig. 23 shows a 
substantial Mach number effect, with incipient vortex angle differences of as 
much as three or four degrees as M c varies from 0.3 to 0.4. Both figures 
indicate that the steady state stall angle based on zero moment slope, 
a ss , is unaffected by sweep angle changes. (However, Figs. 14 through 17 
of Ref. 1 show that the shape of the lift curve changes significantly in the 
direction of lift stall delay with increasing sweep, consistent with Fig. 3 
of Ref. 5 which is based on the data of Ref. 6.) A comparison of the steady 
state stall angle with the zero frequency intercept of the least squares 
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vortex inception line shows a two degree difference at M c = 0.3, and a 
reasonably good correlation at M c = 0.4. These zero frequency and steady 
state values are summarized in Table 3. 


A careful study of Figs. 22 and 23 reveals little or no effect of mean 
angle of attack on the incipient vortex angle. Specifically, at any constant 
value of k c there is no discernable trend of cty with a^; there is only 
scatter. This is an observation that is directly concerned with the effect 
of pitch rate on Oy, which is addressed next. 


A case to be considered in detail is represented by the solid symbols at 
k c = 0.125 in the lower left panel of Fig. 22 for M c = 0.3, a = 10 deg, and 
A = 30 deg. In this instance the vortex inception angle is very nearly 16 
deg for all three values of a^. This is demonstrated in Fig. 24 which 
depicts three motion time histories for = 9, 12, and 15 deg versus 
dimensionless time, t. A horizontal line at Oy = 16 deg represents the 
observed insensitivity of the vortex inception angle to changes in o^, and 
leads directly to a consideration of pitch rate effects. For a sinusoidal 
motion involving instantaneous angle of attack, a, mean angle of attack, 
a M , and amplitude of motion, a as given by Eq . (6), the instantaneous 
dimensional pitch rate is given by 


a = toa cos cot 


<oa 


'sj l-sin^cot 


( 12 ) 


A combination of Eqs. (6) and (12) yields 

a = to ^ - (a-o^)^ (13) 

which can be converted to the dimensionless pitch rate, A, by the formula 
(cf. Ref. 7) 

a = =kc ^ 2 * (a -% )2 (14) 

The values of A for the intersections of the oiy = 16 deg line with the 
three motion curves in Fig. 24 are noted on the figure and are seen to vary 
from 0.0156 to 0.0216. (The minor changes in the actual measured values of 
oty at each k c in Figs. 22 and 23 will result in some scatter but will not 
materially affect the conclusions reached.) 

A direct examination of the effect of mean angle of attack variations on 
oty with reduced frequency effects excluded is afforded by a return to 
Eq. (11). Recalling that this equation is a linear least squares fit, the 
scatter of individual data points may be quantitatively evaluated by 
rewriting the equation in the form 

e “ “v" <Vo~ mk c (15) 
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where e = 0 for an exact linear fit. Values of k c> t q , a v (from Eq. (10)), 
e (from Eq. (15)), and Ay (from Eq. (14)) with a * a y ) are listed in 
Table 4, and e is plotted versus Ay in Fig. 25 for all vortex inception 
events. This figure emphasizes the irregularity of the scatter, and leads 
directly to the observation that the vortex inception angle is independent of 
the pitch rate when k c is fixed. With the exception of a few isolated 
points the error, e, is generally bounded by |e f < 0.5 deg, and never exceeds 
± 1.0 deg. 

A question must be raised at this point on the possible relationship of 
the vortex inception angle, a v , to the various dynamic stall angle concepts 
that have been cited in the literature. As far back as Halfman's 1951 report 
(Ref. 8) and up to the present time, several authors have shown that the 
dynamic stall angle increases with both increasing reduced frequency and 
pitch rate. One of the more complete examinations is that of Gormont (Ref. 
9), based on data and analyses from previous work (Refs. 10, 11). Here it is 
shown that the dynamic moment stall angle, based on 3 C^(dynamic) /3a = 0, 
varies linearly with the square root of the pitch rate, and, for the thick 
V23010-1.58 airfoil examined, approaches the steady state stall angle as the 
zero pitch rate limit is reached. 

To permit a comparison of Gormont 1 s linear fit of the Ref. 10 dynamic 
stall angle results with the current data, the steady-state stall angle was 
removed from the vortex inception angle results of the present experiment 
(a v - a ss , from Tables 3 and 4), and then plotted versus /X in Fig. 26. The 
M c = 0.4 linear fit to Gormont' s thick airfoil moment stall results (lower 
left panel of Fig. 11 in Ref. 9) are shown plotted as dashed lines in the two 
M c = 0.4 panels of Fig. 26. Although there is considerable scatter of the 
data, a linear trend is discernable for a = 10 deg (lower two panels). 

Gormont 's linear fit for moment stall angle starts from the origin and lies 
roughly parallel to, and significantly above, the a y data in the lower 
right panel. If these data were linearly extrapolated to zero pitch rate, 
the intercept value would lie approximately 1.5 deg below the steady state 
stall angle, which is in direct contradiction to the result shown in Table 3 
and in the lower right panel of Fig. 22. This can be resolved by noting that 
a linear behavior of the data is observed for a v - a gs plotted versus /X in 
Fig. 26 or versus A in Fig. 27. In the latter, only the cases for a = 10 deg 
are repeated. The data are now more in keeping with the results in Fig. 22, 
and the original Gormont linear data fit is now a parabola. The apparent 
inconsistency lies in the lack of data at low pitch rates (or reduced 
frequency), and the only clear result from the M c 15 0.4 evidence is that 
leading edge vortex inception precedes moment stall, assuming that both the 
V23010-1.58 airfoil data of Ref. 9 and the current NACA 0012 airfoil data can 
be validly compared in this manner. 
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DISCUSSION OF NORMALIZATION PROCEDURE 


Generally, the purpose of normalization is to consolidate data to a 
common reference. This is done to highlight the nonsimilar aspects of the 
data from case to case and thus to permit a more effective analysis of their 
dependence on the various test parameters. In well behaved small disturbance 
problems, for example, the normalization process coalesces the data to a 
common curve so that they become amenable to a general qualitative analysis. 
In the present investigation, however, the fluctuations in the measured pres- 
sures exceed what is normally considered "small disturbance" levels. In 
realization of this and the fact that much of the data was obtained within 
the stalled flow regime, the following discussion attempts to identify the 
limitations of the normalization procedure used in this report. 

An example of how well the normalization procedure works below stall 
with the current data base is depicted in Fig. 28a in which unsteady data 
obtained from the 30 deg swept wing at a = 8 deg, a M = 0 deg, M c = 0.40, and 
f = 8 Hz (k c = TTcf/M c a 00 = 0.076) are compared with data similarly 
obtained at M c = 0.30, and f = 6 Hz (also at k c = 0.076). This example 
shows that the reduced frequency normalization procedure successfully 
coalesces the linear part (first harmonic component) of the chordwise 
response distribution to a common curve at low aerodynamic load (a M = 0 
deg). Conversely, Fig. 28b shows some deterioration in similarity when 
results for the same combinations of Mach number and frequency are compared 
at = 9 deg. This result indicates that reduced frequency similarity 
deteriorates at high angles of attack. 

The data set reported on herein and in Ref. 1 includes effects due to 
wing sweep angle. Sweep was accounted for in the normalization process by 
utilizing the component of freestream velocity normal to the span as the 
normalizing velocity. This simple approach is the same that has been tra- 
ditionally used in the analysis of steady-state data. An example of the 
effectiveness of this approach to swept wing normalization in similarity 
analysis is shown in Fig. 29a for the zero mean angle of attack case at 
f = 8 Hz, a = 8 deg, and M c = 0.40. Figure 29b shows the same comparison but 
at a M = 9 deg. As in Fig. 28b, it is seen that the higher mean angle of 
attack causes a deteriotat ion in the sweep similarity between normalized 
data . 


On the basis of the examples shown in Figs. 28 and 29, it would appear 
that the present normalization procedure is generally effective in high- 
lighting the nonsimilar aspects of the data, and is thus useful in deter- 
mining the effect of the various test parameters on the oscillatory stall 
phenomenon. 
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It is useful to further examine this normalization scheme by using it to 
make a preliminary evaluation of the effects of f, M c , and A on the com- 
plete unsteady chordwise response of the upper surface at zero mean angle of 
attack. Figures 30 through 32, respectively, show the effect of f, M c , and 
A at the lowest, mean, and highest incidence points of the motion cycle. For 
the range of parameters tested, these figures show that the general shape and 
magnitude of the unsteady chordwise pressure distribution of the zero mean 
angle of attack configuration are in good agreement at the moderate pitching 
amplitude of a = 8 deg. 

The effectiveness of the current normalization procedure in coalescing 
the pressure coefficient distributions helps emphasize the differences in the 
leading edge region results. These differences are especially notable in 
Fig. 31 where the effect of Mach number is shown. Since conclusions about 
the behavior of oscillatory pressures cannot be made on the basis of chord- 
wise pressure distributions alone, individual pressure time histories must 
also be examined to obtain information about the amplitude and phase behavior 
over the entire cycle. An illustration of this is shown in Fig. 33 where 
each column contains cyclically averaged time history traces of the four 
leading edge region pressure coefficients from the same cases that were used 
to produce the chordwise plot comparisons of Figs. 30 through 32. From left 
to right columns, the results in Fig. 33 show the effects of f (corresponding 
to Fig. 30), M (corresponding to Fig. 31) and A (corresponding to Fig. 

32). 


Figure 33 shows that the effect of increasing f on the pressure coeffi- 
cient time history is to increase the lag of the response, while the effect 
of increasing M c is to increase the magnitude of the response at the high 
and low pitching angles of the motion cycle. The rightmost column of Fig. 33 
shows that the effect of A on the pressure coefficient response is minimal 
with the swept configuration slightly leading the unswept case. 

Finally, it is noted that the parameter k c can be varied either by 
changing f or M c . A comparison of the left and middle columns of Fig. 33 
shows that the effect of variation in k c in these figures is dependent on 
whether f or M c is varied. Because of this lack of uniqueness in the 
sensitivity of the response to variations in k c , the parameter analysis in 
this report will focus on the effects of f, M c , and A rather than k c , 

M c , and A . 
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UNSTEADY CHORDWISE PRESSURE DISTRIBUTION ANALYSIS 


Overview 

The chordwise pressure distribution analysis in this section is limited 
to 8 degree pitching amplitude data with special attention to the leading 
edge region. Three examples of the behavior of the chordwise distribution 
are shown in Figs. 34 through 36 for mean angles of attack of 0, 9, and 15 
deg, respectively. Each figure shows the pressure distributions on both 
surfaces at specific points of the motion cycle for both the unswept (left 
column) and swept (right column) wing configurations. The airfoil pitch 
position (relative to mean angle of attack) corresponding to these distribu- 
tions is indicated by the dot on the adjacent waveform between the columns. 
Other test parameters, common to these examples, are M c = 0.30, f = 8 Hz 
(k c = 0.102), and a = 8 deg. 

Before proceeding with the description of the results, it is noted that 
steady state stall occurs at approximately a - 13 deg for M c = 0.30 and at 
approximately a = 10.5 deg for M c = 0.40 (Table 3). Therefore, in terms of 
the present test configurations, dynamic stall penetration occurs at angles 
of attack greater than the mean for a M = 9 deg, either slightly above or 
below the mean for = 12 deg (depending on M c ) and below the mean for 
o^j = 15 deg. This aerodynamic behavior relative to should be borne in 
mind during the discussion of the examples in this section. 

The cases shown in Figs. 34 through 36 are singled out as typical of the 
entire data base for the range of mean angles of attack tested. In addition 
to being a representative zero degree mean angle of attack case. Fig. 34 can 
also be used as a benchmark for evaluating the quality and consistency of the 
data. For example, the upper and lower surface chordwise pressure distribu- 
tions in Fig. 34 appear nearly symmetric with the motion within each column. 
This is consistent with the profile shape of the NACA 0012 airfoil. In par- 
ticular, note how the main features of the upper and lower surface responses 
are reversed at mean angle of attack (a = 0 deg) when the upstroke and down- 
stroke distributions are compared. Similarly note the reversal of the 
surface distributions at the maximum and minimum points of the motion cycle. 
Another test of consistency is to compare the unswept and swept results. As 
expected by the discussion of the previous section, the data appear nearly 
the same throughout the motion cycle. (It is further noted that the unswept 
and swept wing results were generated in two separate tunnel entries.) 

In contrast, the stalled flow penetration case shown in Fig. 35 (o^j = 

9 deg) shows no evidence of symmetry in the response about the mean angle of 
attack. This example also shows that sweep similarity is significantly 
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altered in the higher incidence angle range as the airfoil oscillates in and 
out of stall. In this comparison, the decrease in similarity is especially 
noticeable at peak incidence where the unswept wing has experienced a more 
severe collapse of the suction response than the swept configuration. 

Figure 36 (a M = 15 deg) shows chordwise distributions that were 
generated mostly within the stalled flow region. Here it is seen that the 
differences between the unswept and swept wing results are not as severe, 
giving the appearance that some sweep similarity is regained when the airfoil 
is in deep stall. However, as will be shown in later paragraphs of this 
section, this is not the case since there remains substantial disagreement, 
especially in the forwardmost 2 percent region of the chord, during the 
entire motion cycle. On the basis of Figs. 35 and 36 it appears that the 
dissimilarities due to sweep are strongest when the airfoil load is in 
transition between an unstalled and a stalled response. 


Effect of on Normalization 

The remaining discussion of the chordwise pressure data is focused on 
the sensitivity of the suction (upper) surface distribution to the various 
test parameters with occasional emphasis on the effects of sweep. The cases 
and their behavior highlighted in this section are typical of the entire data 
set . 


In the previous section it was shown that the normalization procedure 
used in the data reduction successfully coalesces the first harmonic compon- 
ent of the zero mean incidence angle cases to a common curve (cf. Figs. 28 
and 29). This coalescence was further demonstrated in Figs. 30 through 32 
for the complete unsteady component of the total response at zero mean angle 
of attack. In these latter three figures comparisons were made showing the 
respective effects of pitching frequency, Mach number, and sweep. 

The impact of mean angle of attack on this generalized behavior is now 
examined by comparing the swept wing results in Figs. 37 (a^ = 0 deg) and 
38 (a M = 9 deg). Each figure, in turn, shows comparisons between two test 
cases in which both the frequency and Mach number are varied simultaneously 
while maintaining a constant reduced frequency of k c - 0.076. The varia- 
tions in frequency and Mach number between cases are the same in each 
figure . 

There is nearly perfect agreement between the two cases over the entire 
cycle at = 0 deg (Fig. 37) confirming the validity of reduced frequency 
scaling at low load. However, when a similar comparison is made at = 

9 deg (Fig. 38), the two cases are in disagreement. This deterioration in 
generalized behavior occurs over the entire high angle of attack range of the 
motion cycle with maximum disagreement at peak a where the unsteady separ- 
ation effects spread over the entire chord. Conversely, it is noted that 
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similarity in the chordwise response is regained at the minimum point of the 
cycle where a * 1.0 deg (Fig. 38). The disagreement in the chordwise 
responses over most of the cycle is due to the difference in Mach number 
between the two cases as illustrated in the following discussion. 

The sensitivity of the data at a M - 9 deg to variations in the reduced 

frequency is examined in the next two figures for both the unswept (left 

column) and swept (right column) wings at the mean, maximum, and minimum 

points of the motion cycle. The objective here is to show the influence of 

k c on the data by separately varying the pitching frequency (Fig. 39) and 

Mach number (Fig. 40). In Fig. 39, the variation in k £ (due to f) is three 

times that of Fig. 40 (due to M Q ) and yet the degree and trend of the 

change in the chordwise response between cases in each figure is about the 

same for both the unswept and swept configurations. These figures show that 

the trend of the pressure response distribution, although different in 

magnitude, is generally consistent with the variation in k £ (due to changes 

in either f or M„) and verify the dominant influence of M„ on k„ . 

c J c c 

A further inspection of Figs. 39 and 40 also shows that although the 
influence of k c (due to changes in either f or M c ) is generally the same 
between unswept and swept results, there are regions of the motion cycle 
where k c has the opposite effect. (For example, near peak angle of attack 
in Fig. 40, an increase in k c corresponds to a decrease in suction for the 
unswept response and an increase in suction for the swept response in the 
forward 10 percent region of the chord.) 

This difference in the effect of k c on the response at specific points 
of the cycle is partly due to the difference in the phasing of the unsteady 
load response relative to the pitching motion between the unswept and swept 
configurations. As a result, the effect of k c shows up differently at 
various points of the cycle, especially in the transition region (between 
stalled and unstalled flow) where the airfoil surface pressure response 
appears to be most sensitive to parameter variations. Examples illustrating 
the sensitivity of the phasing of the leading edge region pressure response 
are presented next. 


Leading Edge Region Time Histories 

An important aspect in the parameter ic evaluation of chordwise load 
behavior is the analysis of the pressure time histories at various points 
along the chord. The objective here is to show representative examples from 
the data base and to highlight some of the main features in the behavior of 
the pressure time histories to variations in the test parameters. The time 
histories shown in Figs. 41 through 43 each contain three-harmonic component 
reconstructions of the leading edge region responses at the 0.4, 1.0, and 1.9 
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percent chord positions of the suction surface. (The choice of three-har- 
monic component plots rather than the total response is to keep the para- 
metric analysis relatively simple while not giving up any essential 
features.) Each figure is arranged in 3 columns corresponding to the 0, 9, 
and 15 degree mean incidence angle configurations. Each column is comprised 
of the blade pitching time history at the top followed by the three forward- 
most suction surface time histories at x = 0.004, 0.010, and 0.019. 

Figure 41 shows the effect of frequency at M c = 0.30, A = 0 deg, and 
a = 8 deg. In this figure, it is seen that a frequency increase from 4 
to 8 Hz (or k c increasing from 0.051 to 0.102) has only a minor impact, 
showing up mainly as a slight increase in the phase lag (right shift) of the 

responses at = 0 deg and 9 deg. At = 15 deg, however, this 

increase in frequency causes a substantial increase in the phase lag of the 
response during the upstroke as well as a significant dimunition in the 
response level during the downstroke. Figure 44 is a companion to Fig. 41 
where the 9 deg (left column) and 15 deg (right column) mean incidence angle 
time histories from the leading edge region are repeated with additional time 
histories further down along the chord at x = 0.045, 0.073, 0.114, and 0.268. 
This figure shows little variation in the relative phasing throughout the 

forwardmost quarterchord region. This is most easily seen during the first 

quarter cycle of the - 15 deg configuration. 

Figures 42 and 45 show the effect of Mach number at f = 8 Hz, A = 0 deg, 
and a - 8 deg using the same comparison approach as in Figs. 41 and 44. Here 
Fig. 42 shows that an increase in Mach number has a greater effect on the 
leading edge region response than does the variation in the frequency (cf. 
Fig. 41). This sensitivity of the leading edge region response to the Mach 
number is especially significant at = 9 deg and to a lesser extent at 
a M = 15 deg. In particular, it is noted that the increase in Mach number 
lowers the stall penetration angle and reduces the sharpness of the dynamic 
stall event. Figure 45 shows that this sensitivity of the aerodynamic 
response to M c persists at least to x = 0.268. 

Examples of the effect of sweep on the time history behavior of the 
chordwise response are shown next in Figs. 43 and 46 at f = 8 Hz, M c = 

0.30, and a - 8 deg. The zero degree mean angle of attack results in Fig. 43 
show that sweep effects are essentially negligible throughout the cycle all 
way to the forwardmost station at x = 0.004. Conversely, sweep signficantly 
alters the leading edge region response when = 9 and 15 deg. In the 
9 deg mean angle of attack case, sweep effects remain fairly significant 
through the quarterchord region while the 15 deg mean angle of attack case 
shows a substantial decline in the influence of sweep beyond x ~ 0.045. The 
most notable effect of sweep on these results is in the increased phase lag 
of the load and the associated delay in the stall event. This is especially 
clear when “M " 9 deg. In particular, note how the peak response during 
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the upstroke is delayed to later in the cycle for A = 30 deg (dashed line) 
relative to A = 0 deg (solid line). 


Parameter Analysis of Chordwise Pressures 

The effect of mean angle of attack on the suction surface pressure 
response and on the influence of sweep at specific angles of attack is 
examined next. As examples. Figs. 47 and 48 show this effect at a = 7.4 deg 
and 16.6 deg in the upstroke portion of the cycle. Each figure shows results 
at M c = 0.30 (left column) and at M £ = 0.40 (right column) along with 
accompanying simple waveforms indicating the point in the motion cycle where 
the chordwise distributions were taken. Figure 47 contains results taken 
from the portion of the cycle where the unsteady aerodynamic activity is not 
in transition and is within the potential flow range of the motion. It is 
seen that the unswept and swept wing chordwise distributions in this part of 
the cycle are each nearly similar over the range of mean angle of attack 
shown at Mach numbers of 0.30 and 0.40. This behavior, although not 
completely generalized, is encouraging from a data utilization standpoint in 
that this response distribution appears applicable (at this point of the 
motion cycle) over the entire range of mean incidence angles between = 

9 deg and = 15 deg. With the exception of the upper left distributions, 
Fig. 48 shows that this similarity of the unswept and swept wing responses 
with respect to mean angle of attack tends to persist into the separated flow 
region. In particular, the chordwise distributions for each Mach number in 

Fig. 48 are very similar at ct^ = 12 deg and 15 deg while noticeably differ- 
ent at = 9 deg where the angle of attack value used in this comparison 

is near its peak. (It is again pointed out here that the jaggedness in the 

leading edge region response is an anomaly of the data which was explained 
earlier.) It is also noted in this figure that the effects of Mach number 
and sweep are substantial at this high value of angle of attack (16.6 deg), 
while the effect of Mach number at a = 7.4 deg (cf. Fig. 47) is nearly 
imperceptible. 

A further analysis of the M c = 0.30 results corresponding to the 
previous two figures is presented in Figs. 49 and 50. As in Figs. 47 and 48, 
these figures show the effects of sweep and mean angle of attack on the 
chordwise distribution at various points of the motion cycle. Figure 49 
shows the chordwise distributions at instantaneous angles of attack of a = 9, 
11, 13, and 15 deg during upstroke (in columns from left to right) thus fill- 
ing the gap between the results of Figs. 47 and 48 at 7.4 deg and 16.6 deg 
incidence, respectively. It is seen that the effect of varying on the 
shape of the distribution (within each column) is relatively small over this 
range of angle of attack (within each row) for both the unswept and swept 
cases, implying that significant stall penetration has not yet taken place. 

It is also noted that the effect of sweep in reducing suction is consistent 
throughout this portion of the motion cycle. 


26 



A similar analysis is presented in Fig. 50 for the downstroke portion of 
the cycle over the same angle of attack range (where the plots are arranged 
following the motion from left to right). In this region the effects of 
sweep and mean angle of attack are each more significant than in Fig. 49. 

This is especially evident at a = 15 deg (extreme left column) where flow 
reattachment is still underway. The sensitivity of the chordwise response to 
changes in mean angle of attack during the downs troke portion of the motion 
cycle occurs mainly because the aerodynamic activity is in varying stages of 
recovery from separation. This variation in the level of aerodynamic activ- 
ity depends both on how deeply the stalled flow region was penetrated by the 
airfoil and on where the airfoil is in the motion cycle relative to peak 
incidence. 


Effect of Sweep on Peak Angle of Attack Response 

The effect of sweep on the peak angle of attack region response is now 
examined for a M = 9, 12 and 15 deg. Figures 51 and 52 show the sensitivity 
of the chordwise response to sweep in 0.5 degree increments within 1.5 deg of 
peak angle for each mean angle of attack case. Figure 51 shows the upstroke 
portion of the study while Fig. 52 shows the downstroke portion. 

The effect of sweep on the 9 degree mean angle of attack configuration 
is examined first (see top row of Figs. 51 and 52). At a = 15.5 deg in the 
top row in Fig. 51, the effect of sweep is still moderate as the airfoil 
motion enters the stall transition region of the cycle. A half degree later 
at a = 16.0 deg, there is beginning evidence of collapse in suction for the 
unswept wing. This collapse progresses rapidly to maximum severity at peak 
angle while the loss in suction for the swept wing occurs in a more orderly 
fashion and does not reach the same intensity. As a result of this differ- 
ence in the rate and magnitude of suction collapse, the unswept and swept 
wing chordwise load distributions are signficantly different at peak angle of 
attack. Going on to the top row in Fig. 52 where the downstroke portion is 
presented, the change in both the unswept and swept wing responses becomes 
less severe with the change in angle of attack and the chordwise distribu- 
tions appear to once again coalesce as a reaches 15.5 degrees. 

The results for the 12 degree mean angle of attack configuration are 
shown in the middle row of both figures. In this case the motion is well 
into the stall region throughout the range of angles shown and the pressure 
distribution is still exhibiting the after-effects of transition between 
unstalled and- stalled flow. This behavior is made more apparent by the 
extensive disagreement between the unswept and swept airfoil pressure distri- 
butions throughout the upstroke portion of the motion near peak angle of 
attack (Fig. 51). This sweep-induced difference becomes less severe as the 
motion passes through the peak and into the downstroke region (Fig. 52). 
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However, the strong disagreement at the leading edge persists throughout this 
range of angles. As shown earlier in the discussion of response time 
histories, this behavior is partly a consequence of the difference in the 
phasing of the recovery process between the unswept and swept wing responses 
relative to the airfoil pitching motion. 

The behavior of the 15 degree mean angle of attack configuration is 
shown in the bottom row of Figs. 51 and 52. In this case, the chordwise 
distributions all occur very deep in the stall region, where flow separation 
is fully achieved, and are thus relatively stable to small changes in angle. 
At a = 21.5 deg in Fig. 51, it is seen that the transition into fully 
separated flow is nearly completed. Once the transition is fully achieved, 
the unswept and swept distributions appear to coalesce and as the motion 
passes through peak angle of attack, the load distributions appear less 
dependent on sweep. It should be noted, however, that this latter 
observation can be misleading if it is not combined with a study of the 
corresponding time histories of Fig. 43 where it is seen that the close 
agreement between unswept and swept results is only transitory during the 
downstroke at x = 0.33 (right column). 
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DISCUSSION OF INTEGRATED LOAD RESULTS 


The integrated load results (lift, drag, and moment) reported in Ref. 1 
were computed from the chordwise response data of this report. Hysteresis 
loop comparisons were made therein to show the effect of sweep on the 
unsteady lift, drag, and moment components of the total force response. A 
particular feature that remained unexplained in the results is the substan- 
tial displacement in the lift curve below stall between the unswept and swept 
wing results, especially at high angles of attack. A perspective on this 
displacement can be obtained from Figs. 53 and 54. Figure 53 shows that the 
normalized steady-state lift curves of the unswept and swept NACA 0012 air- 
foils are nearly coalesced below stall. This is not the case' when viewing 
the unsteady comparisons in Fig. 54. Here, it is seen that the upstroke 
portions of the unswept lift responses below stall are significantly and 
increasingly displaced from one another as a M increases. (Figure 54 is a 
reproduction of results from Fig. 34 of Ref. 1 showing the effect of sweep 
on the lift response for four mean angles of attack tested at M c = 0.40, 
f = Hz, and a = 8 deg.) 

It has been found that this displacement of the lift reponse is mainly 
associated with the first harmonic component. This was discovered serendipi- 
tously while attempting to superimpose the unsteady parts of the unswept and 
swept wing lift responses. In order to achieve this superposition effec- 
tively, the time mean components of the lift responses were removed in an 
attempt to bring the loops to closer alignment than shown in Fig. 54. The 
result of doing this is shown in the top row of Fig. 55 for the 15 degree 
mean angle of attack case where plot 55a is repeated from Fig. 54 and plot 
55b is stripped of the time mean component. The resultant lack of change in 
the relative positioning of these loops indicates that, for this particular 
comparison, the entire displacement is a sweep-induced unsteady phenomenon. 

A harmonic analysis of these loops was performed to determine the aspect of 
the data that accounts for this displacement. It was found that the ultimate 
shape, displacement, and rotation of the lift loops of Fig. 55a are primarily 
governed by the first harmonic component ellipse. This is clearly illus- 
trated by the first harmonic component comparison in plot 55c where the mid- 
portion of the AC l upstroke displacement of the complete response (plot 
55a) is accounted for. This first harmonic result shows that the primary 
influence of sweep on the total lift response is through the phasing between 
the first harmonic component of the lift response and the blade motion, and 
through the ratio of the first harmonic component amplitude of the response 
to the airfoil pitching amplitude. The remaining comparisons in Fig. 55 
(plots d through h) show how the addition of higher harmonic terms contrib- 
utes in evolving the complete unsteady response without significant influence 
on the AC^ displacement at the midpoint of the upstroke. From this 
example, it is seen that the main importance of the higher harmonic terms is 
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in shaping the loops rather than affecting their orientation. A comparison 
of plot 55d and plot 55b shows that most of the remaining features of the 
loops are due to the second harmonic term. 

The analysis of the moment responses for the same unswept and swept 
cases shown in Fig. 55 is presented in Fig. 56. Unlike the lift result, it 
is seen that the removal of the nonoscillatory terms from the total values 
tends to bring the unswept and swept responses to closer proximity but with 
some displacement remaining. As in the case of the lift result the subse- 
quent harmonic analysis of the moment response shows that the remaining dis- 
placement is due to the influence of sweep on the unsteady portion of the 
response and that the ultimate orientations are again largely determined by 
the first harmonic ellipses. In this case the influence of the second 
harmonic term (plot 56d) on the final shape is clearly evident. 

Several additional sweep comparisons were analyzed and variations were 
observed in the source (unsteady or time mean) of the loop displacement from 
case to case. An example that illustrates this variation in the sweep influ- 
ence (relative to that shown in Figs. 55 and 56) is presented next in Figs. 

57 and 58. 

The unswept and swept cases shown in these figures are the same in every 
respect to Figs. 55 and 56, except that - 9 deg. The results of the 
harmonic analysis of the lift in Fig. 57 is basically the same as before; 
however. Fig. 58 shows that, in this case, the displacement of the moment 
response is entirely accounted for by the nonoscillatory terras. As a result, 
note how the first harmonic ellipses are nearly coincident. 

Taken together, the harmonic analysis presented in Figs. 55 through 58 
show that the loop displacements range from being entirely dependent on the 
unsteady influence (cf. Fig. 55) to entirely dependent on the time mean 
influence (cf. Fig. 58). Therefore, what once appeared to be anomalous 
behavior (because of the inconsistent result obtained by removal of the time 
mean component from the total response) has been shown to be amenable to 
analysis . 
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SUMMARY OF RESULTS 


An analysis was performed on unsteady aerodynamic pressure data from a 
NACA 0012 tunnel-spanning wing oscillating in the unswept and swept configu- 
ration. The range of parameters tested were: sweep, 0 and 30 degrees; 

pitching amplitude, 8 and 10 degrees; mean angle of attack, 0 to 15 degrees; 
Mach number normal to the wing span, 0.30 and 0.40; pitching frequency, 4 to 
11 Hertz. The corresponding reduced frequency range of the data was .038 to 
.125. 


The main subject areas of this report are: the analysis of the chord- 

wise wave propagation speed, the analysis of the upper surface chordwise 
pressure distribution, and a brief harmonic analysis of the integrated lift 
response. The results of this study are summarized below. Items with 
asterisks are regarded as primary findings. 


Pressure Wave Analysis 

• Mean angle of attack has little or no effect on wave speed, which in 
all cases increases uniformly with reduced frequency (Figs. 16 through 
20). Motion amplitude also has little or no effect on wave speed (Fig. 
16). Mach number has a small but consistent effect on wave speed, with 
higher Mach numbers yielding slightly smaller wave speeds (Fig. 17). 

* • Sweep angle has a dominant effect on wave speed. The overall wave 

speed for the unswept wing is consistently greater than that for the 
swept wing by a factor of 1.5 to 2.0 (Figs. 18 et seq.). This 
represents a major failure above stall for the cosine law normalization 
which has been shown to be consistently valid below stall. 

• Previous observations of the effect of increasing sweep to delay stall 
are confirmed in the present work. In addition, the vortex inception 
angle increases slowly with k c (Fig. 22). 

* • Mach number has a dominant effect on vortex inception angle and static 

stall angle, with a decrease in both angles as M increases (Fig. 

23). 

• Vortex inception is substantially independent of amplitude of motion 
and of mean angle of attack (Figs. 22, 23). 

• Local wave speed differences associated with sweep are confined to the 
region of the blade aft of the 15 percent chord. Wave speeds are 
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substantially the same for both sweep angles forward of this chord loca- 
tion (Figs. 19, 20). The variation of local wave speed with chordwise 
position shows an initial rapid increase near the leading edge followed 
by a leveling off of wave speed over the aft region of the chord (Fig. 
21). Overall wave speed (averaged over the first 15 percent of the 
chord) agrees well with previous results from Ref. 4 (Fig. 15). 

* • The vortex inception angle is independent of pitch rate when reduced 

frequency is fixed (Fig. 25). 

Chordwise Pressures Analysis 

• Sweep effects on the upper surface pressure distribution are more 
prevalent when the aerodynamic response is in transition between 
unstalled and stalled conditions. 

• The sensitivity of .the pressure responses to sweep is also strongest 
when the wing angle of attack is near peak value for mean angles of 
attack of 9, 12, and 15 degrees (Figs. 34, 35, and 36). 

* • Sweep similarity by means of the cosine law (which prevails over the 

angle of attack range of ± 8 degrees about zero mean angle) deteriorates 
when the mean angle of attack reaches 9 degrees and the motion pene- 
trates stall (Figs. 43 and 46). 

• The effects of sweep very near peak angle appear stronger on the 
approach to than on recession from peak angle when = 9 deg while 
the opposite appears to be the case when = 12 deg. The effects of 
sweep at a^. = 15 deg, on the other hand, are relatively small near 
peak angle (Figs. 51 and 52). 

* • The reduced frequency parameter is not a reliable indicator of the 

unsteady aerodynamic condition at high angles of attack. This is 
because the same numerical value of reduced frequency can be obtained 
for two different combinations of M c and f in which the M c effect on 
the flow dominates the effect of f (Figs. 37 and 38). 

Integrated Load Analysis 

* • The main influence of sweep on the total integrated lift and moment 

responses is very strongly determined by the first harmonic term. This 
result implies that the effect of sweep on the unsteady response is 
amenable to analysis (Fig. 55). 
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APPENDIX 


DISCUSSION OF MEASUREMENT ERROR 


A minor measurement error was detected during the present analysis that 
occurs throughout the data set. This error neither affects the discussion 
presented herein nor does it alter the conclusions previously reported in 
Ref. 1 (based on the same data set). 

The main evidence of the error is the jaggedness that occurs in both the 
upper surface chordwise pressure distribution plots and the constant pressure 
contours of the chordwise wave speed analysis. For the chordwise distribu- 
tion plots, this anomaly ranges from slight to relatively large depending on 
the pressure loading condition. In the case of the constant pressure con- 
tours the effect is generally minor and, at worst, is more a nuisance than a 
detriment to the evaluation of the wave speed. 

A close inspection of the first harmonic component phase angle distribu- 
tion along the chord in Fig. 4 shows that the jaggedness (seen in Fig. 5 for 
example) is principally related to a phasing error between consecutive pres- 
sure responses along the chord. In particular, it is noted that the phase 
angle plots near the leading edge appear to be scattered about a mean trend. 
This behavior of the data set was found to be caused by a misalignment of the 
tape heads of the analog tape recorders. Both the read and write functions 
were configured to utilize two heads during record and playback where the odd 
and even numbered channels (sequentially assigned along the chord) were 
segregated to separate heads. 

Under these circumstances, the time lag between odd and even channels 
due to head misalignment caused an odd-even channel phase error. This phase 
error, A , is related to the oscillation period, T, and tape speed, S, by the 
equat ion 

-SiL . Ax = At. fAt = fAX 
360* T S 

where Ay represents the error displacement between the heads. A typical 
value of Ay equal to .003 in. produces a phase angle error of 5 deg in the 
first harmonic at f = 8 Hz . 

The entire data set was uniformly corrected for the computations made in 
the wave speed analysis, but the corrections were confined to that section 
only. All other figures presented in this report and the data tables in Ref. 
3 remain uncorrected to preserve consistency with the results of Refs. 1 
and 2 . 
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TABLE 1 


UNSTEADY TEST MATRIX FOR NACA0012 AIRFOIL 


Airfoil 

A 

(deg) 

a 

(deg) 

a M 

(deg) 

M. 

M c 

f(cps) 

4 

5.33 

6 

8 

10 

10.67 

NACAOOl 2 

0 


S 

0 

.3 

.3 

X 



X 

X 








ii 

.4 

.4 

X 



X 

X 







Jj| 

9 

.3 

.3 

X 



X 

X 







B 

it 

.4 

.4 

X 



X 

X 







■ 

12 

.3 

.3 

X 



X 

X 







SB 

II 

.4 

.4 

X 



X 

X 








15 

.3 

.3 

X 



X 

X 








II 

.4 

.4 

X 



X 

X 






Ii 

I 

0 

.3 

.3 

X 



X 

X 






mBg? 


II 

.4 

.4 

X 



X 

X 








9 

.3 

.3 

X 



X 

X 









.4 

.4 

X 



X 

X 







■ 


.3 

.3 

X 



X 

X 






■ 

| 


.4 

.4 

X 



X 

X 


NACAOOl 2 

30 

8 

0 

.346 

.3 

X 

■ 

X 

X 

X 








n 

.462 

.4 

X 

1 


X 

X 

X 







9 

.346 

.3 

X 

II 

X 

X 

X 








II 

.462 

.4 

X 

1 


X 

X 

X 







12 

.346 

.3 

X 

H 

X 

X 

X 








II 

.462 

.4 

X 

m 


X 

X 

X 







15 

.346 

.3 

X 

M 

■a 

X 

X 








II 

.462 

.4 

X 

mm 

n 

X 

X 

X 





10 

9 

.346 

.3 

X 


n 

X 

X 








n 

.462 

.4 

X 

X 

■ 

X 

X 

X 







12 

.346 

.3 

X 


■I 

X 

X 








n 

.462 

.4 

X 

mm 


■3 

X 

X 







15 

.346 

.3 

X 

H 

n 

u 

X 








ii 

.462 

.4 

X 

D 

H 

n 

X 

X 
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Notes : 

1. All values displayed were run at least once. 

2. Horizontal arrows on right column denote nominal values 
used for unswept runs. 

3. Boxed numbers denote nominal values used for swept runs. 

4. Diagonal arrows indicate matched values of k c . 
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TABLE 3 




Least Squares 
intercept at 

V° 


Steady State 
Based on 
dC M /da=0 

M c 

a 


A=0° 

A=30 ° 


A=0° 

A=30 ° 

.3 

8° 


14.4° 

15.0° 


12 .9° 

12 .8° 

10° 


14.5° 

14.5° 

.4 

8° 


10.3° 

11 .4° 


10.4° 

10.4° 

10° 


9.9° 

►— * 
o 

On 

o 
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TABLE 4 




k 

T 

a 

£ 

A 


, M . 
(deg) 

c 

O 

(deg) 

(deg) 

V 


9 

.0498 

.137 

15.07 

.423 

.0045 

A=0 deg 


.0993 

.125 

14.66 

-.202 

.0098 



.1245 

.123 

14.58 

-.391 

.0125 

M =0.3 

12 

.0497 

.052 

14.57 

-.077 

.0066 



.1001 

.057 

14.80 

-.065 

.0131 

ff=8 deg 


.1249 

.066 

15.22 

.247 

.0160 


15 

.0497 

.987 

14.35 

-.297 

.0069 



.1002 

.000 

15.00 

.134 

.0140 



.1248 

.004 

15.20 

.228 

.0174 


9 

.0499 

.094 

14.57 

-.334 

.0072 

A=0 deg 


.0997 

.104 

15.08 

-.195 

.0138 



.1250 

.110 

15.37 

-.093 

.0168 

M c =3.0 

12 

.0497 

.051 

15.15 

.248 

.0082 



.1006 

.061 

15.74 

.459 

.0163 

ff=10 deg 


.1265 

.055 

15.39 

-.085 

.0208 


9 

.0374 

.050 

11 .47 

.908 

.0050 

A=0 deg 


.0742 

.033 

10.65 

-.213 

.0101 



.0929 

.045 

11 .23 

.214 

.0125 

M c =0.4 

12 

.0374 

.969 

10.45 

-.112 

.0051 



.0749 

.986 

11.30 

.431 

.0104 

5=8 deg 


.0939 

.991 

11 .55 

.526 

.0131 


15 

.0373 

.889 

9.86 

-.701 

.0040 



.0751 

.902 

10.38 

-.491 

.0086 



.0935 

.904 

10.46 

-.561 

.0108 


9 

.0374 

.026 

10.63 

-.054 

.0064 

A=0 deg 


.0746 

.033 

11.06 

-.369 

.0127 



.0934 

.042 

11 .61 

-.196 

.0157 

M c =0 .4 

12 

.0376 

.979 

10.68 

-.008 

.0065 



.0749 

.000 

12.00 

.565 

.0131 

5=10 deg 


.0935 

.998 

11 .87 

.062 

.0163 


9 

.0502 

.138 

15.10 

-.185 

.0045 

A=30 deg 


.0756 

.144 

15.29 

-.148 

.0065 



.1003 

.149 

15.44 

-.147 

.0083 

M c =0.3 


.1259 

.145 

15.32 

-.421 

.0108 


12 

.0498 

.065 

15.18 

-.103 

.0064 

5=8 deg 


.0656 

.074 

15.59 

.212 

.0082 



.0750 

.071 

15.45 

.015 

.0095 



.0996 

.079 

15.81 

.227 

.0122 



.1244 

.074 

15.59 

-.142 

.0155 



.0503 

.060 

14.94 

-.346 

.0065 



.1003 

.077 

15.72 

.133 

.0124 


15 

.0503 

.013 

15.65 

.364 

.0070 



.0758 

.009 

15.45 

.011 

.0106 



.1000 

.018 

15.90 

.315 

.0139 



.1246 

.019 

15.95 

.216 

.0173 


a vo Bl4 * 43# 

M c =4.330° 


V 14 * 53 ' 

M c =7 .454° 


%o =10 - 260 

M c =8.177° 


a =9.94° 
vo 

M c =20.023° 


a vo~ 14 ,9 ®° 
M c =6.026° 
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TABLE 4 (Concluded) 



, a M 
( deg) 

k c 

T 

O 

(deg) 

e 

(deg) 

A v 


9 

.0498 

.102 

14.98 

-.111 

.0070 

A=30 deg 


.0750 

.110 

15.37 

-.013 

.0101 



.0994 

.114 

15.57 

-.095 

.0131 

M =0.3 


.1241 

.115 

15.61 

-.341 

.0163 


12 

.0497 

.049 

15.03 

-.060 

.0083 

3—10 deg 


.0754 

.050 

15.09 

-.297 

.0125 



.0995 

.060 

15.68 

.014 

.0162 



.1241 

.066 

16 .03 

.079 

.0198 


15 

.0499 

.005 

15.31 

.218 

.0087 



.0754 

.008 

15.50 

.113 

.0131 



.0995 

.010 

15.63 

-.036 

.0173 



.1239 

.018 

16.13 

.181 

.0215 



.0983 

.016 

16.00 

.348 

.0171 


9 

.0377 

.051 

11 .52 

-.269 

.0050 

A=30 deg 


.0502 

.055 

11 .71 

-.223 

.0066 



.0754 

.050 

11 .47 

-.752 

.0100 

M =0 .4 


.0945 

.069 

12.36 

.033 

.0107 



.1002 

.068 

12 .32 

-.187 

.0127 

3=8 deg 

12 

.0376 

.003 

12 .15 

.362 

.0052 



.0496 

.010 

12.50 

.574 

.0069 



.0745 

.014 

12.70 

.488 

.0104 



.0933 

.016 

12.80 

.372 

.0130 



.0992 

.015 

12 .75 

.254 

.0138 



.0377 

.998 

11 .90 

.111 

.0053 



.0753 

.002 

12.10 

-.121 

.0105 


15 

.0376 

.929 

11.55 

-.238 

.0047 



.0499 

.934 

11 .78 

-.149 

.0064 



.0751 

.941 

12.10 

-.119 

.0098 



.0944 

.954 

12 .72 

.280 

.0126 



.0997 

.943 

12.20 

-.301 

.0130 


9 

.0371 

.035 

11 .18 

-.107 

.0063 

A=30 deg 


.0495 

.042 

11.61 

.096 

.0083 



.0743 

.045 

11 .79 

-.176 

.0125 

M =0 .4 


.0929 

.050 

12.09 

-.216 

.0154 



.0990 

.052 

12.21 

-.207 

.0164 

3=10 deg 

12 

.0375 

.989 

11 .31 

.016 

.0065 



.0497 

.992 

11.50 

-.017 

.0087 



.0746 

.003 

12.19 

.218 

.0130 



.0935 

.000 

12.00 

-.317 

.0163 



.0992 

.010 

12.63 

.209 

.0173 


15 

.0374 

.935 

11 .03 

-.263 

.0060 



.0494 

.949 

11 .85 

.338 

.0082 



.0740 

.953 

12.09 

.129 

.0124 



.0934 

.957 

12.33 

.015 

.0157 



.0991 

.963 

12.70 

.281 

.0168 


V 14 ' 51 ' 

M =11 .576" 
c 


%o =11 - 36 “ 
M c =ll .487° 


a vo =l°.61" 

M =18 .254° 
c 
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Fig. 1 Airfoil Cross Section Showing Chordwise Measuring Stations Along Reference Line 



CHORDWISE COMPONENT OF 
FREESTREAM VELOCITY 


C 0 - C N sin fl - Cc cos a 
C L • C N COS a * Cq SIN a 


Fig. 2 Airfoil Force Schematic 


41 


UNSTEADY PRESSURE COEFFICIENT C p (*) UNSTEADY PRESSURE COEFFICIENT C p (*) 
(SUCTION SURFACE) (SUCTION SURFACE) ^ 


- 




MEAN INCIDENCE 

{ 




M c = 0.30 

— 




M c = 0.40 

i L 

1 

1 

1 

n n n n 1 



DIMENSIONLESS CHORDWISE POSITION, x 


Fig. 3 Example of Chordwise Pressure Distribution Comparison — Continuous Curves 
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a M = 9deg 



DIMENSIONLESS CHORDWISE POSITION, * 


Fig. 4 Example of Chordwlse Pressure Distribution Comparison — Discrete Points 
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UNSTEADY PRESSURE COEFFICIENT, C 



DIMENSIONLESS CHORDWISE POSITION, x 


Fig. 5 Example of Chordwlse Pressure Distribution Format for Comparison 

of Adjacent Columns 
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Fig. 6 Example of Chordwise Pressure Distribution Format for Trend Analysis 


THREE HARMONIC COMPONENT RECONSTRUCTION OF 
PRESSURE COEFFICIENT. C D (x) 


M c = 0.30 M c = 0.40 

a M = 9 deg a M = 15deg 

ALPHA 

4 

X = 0.004 o 
-4 

4 

X = 0.0t0 0 

-4 








Fig. 7 Example of Pressure Time History Format for Comparison or Trend Analysis 
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Pressure Carpet Plot, Leading Edge Region 
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Ml .167 .334 .Ml .669 .931 1.M9 

0 I BENS I0NLE8S TIRE. FRACTION OF PERIOD , T 



Fig. 10 Pressure Carpet Plot, Entire Chord 
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Fig. 12 Pressure Contour Map with Wave Speed Definition, Near Leading Edge 
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Fig. 13 Illustration of Segmented Wave Speed Near Leading Edge 
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Fig. 15 Variation of Normalized Wave Speeds with Frequency. _ 
(In Ref. 4 V Wq2 was Denoted V w<| and V Wq3 was Denoted V w ) 
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Fig. 17 Effect of Mach Number on Overall Wave Speed 
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Fig. 18 Effect of Sweep Angle on Overall Wave Speed 
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Fig. 21 Chordwise Variation of Local Wave Speed for =15 deg, M 
and a = 8 deg 
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Fig. 23 Effect of Mach Number on Vortex Inception Angle Based on Initial Wave 
Propagation 
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Fig. 27 Correlation of Relative Vortex Inception Angle with A v 
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DIMENSIONLESS CHORDWISE POSITION, x 


Fig. 30 Effect of Frequency on Complete Unsteady Chordwlse Pressure Distribution at the 
Lowest v Mean, and Highest Angle of Attack Points of Motion Cycte 
(c*m * o deg, M c = 0.30, 5 = 8 deg, A = 0 deg) 
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Fig. 31 Effect of Mach Number on Complete Unsteady Chordwfse Pressure Distribution at 
the Lowest, Mean and Highest Angle of Attack Points of Motion Cycle 
(c*M -0 deg, f =8 Hz, a = 8 deg, A = 0 deg) 
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Fig. 32 Effect of Sweep on Complete Unsteady Chordwise Pressure Distribution at the 
Lowest, Mean, and Highest Angle of Attack Points of Motion Cycle 
=0 deg, k c =0.100 (f = 8 Hz, M c =0.30), end a = 8 deg 
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Fig. 33 Parameter Sensitivity of Upper Surface, Leading Edge Region Pressure Coefficient Time Histories 

at = 0 deg and « = 8 deg 
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Fig. 34 Chordwise Pressure Results for the Unswept and Swept Configurations at 
=0 deg, k c =0.100 (f =8 Hz, M c =0.30), and a = 8 deg 
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Fig. 35 Chordwlse Pressure Results for the Unswept and Swept Configurations at = 9 deg, 

k c =0.100 (f =8 Hz, M c = 0.30), and a =8 deg 
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Fig. 36 Chordwise Pressure Results for the Unswept and Swept Configurations at =15 deg, 

k c = 0.100 (f = 8 Hz, M c = 0.30), and a = 8 deg 
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Fig. 38 Upper Surface Chordwise Pressure Behavior of Swept Airfoil at =9 deg, 

k c = 0.075, and a =8 deg 
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Fig. 39 Effect of Frequency on Upper Surface Response at Mean, Maximum, and Minimum 
Points of Motion Cycle =9 deg, M c = 0.30, 5 = 8 deg) 


78 



DIMENSIONLESS CHORDWISE POSITION, X 


Fig. 40 Effect of Mach Number on Upper Surface Response at Mean, Maximum, and 
Minimum Points of Motion Cycle ( a m =9 deg,f =8 Hz, 5=8 deg) 
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Fig. 41 Effect 
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Fig. 42 Effect of Mach Number on Upper Surface Leading Edge Region Time Histories (f = 8 Hz, a = 8 deg, A = 0 deg) 
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Fig. 43 Effect of Sweep on Upper Surface Leading Edge Region Time Histories 
k c = 0.100 (f = 8 Hz, M c = 0.30) and Ti = 8 deg 
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Fig. 45 Effect of Mach Number on the 
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Fig. 47 Effect of on Sweep Dependence of Upper Surface Response at a = 7.4 deg 
Upstroke (f =8 Hz, a = 6 deg) 
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Fig. 48 Effect of ckm on Sweep of Dependence of Upper Surface Response at a = 16.6 deg 
Upstroke (f = 8 Hz, a = 8 deg) 
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Fig. 53 Steady-State Lift Coefficient for NACA 0012 Airfoil at M c = 0.40 
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Fig. 55 Harmonic Analysis of the Sweep Effect on the Lift Response at a M =15 deg, 
k c = 0.075 (f = 8 Hz, M c = 0.40), and a =8 deg 
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Fig. 56 Harmonic Analysis of the Sweep Effect on the Moment Response at = 15 deg, 

k c =0.075 (f =8 Hz, M c =0.40), and a =8 deg 
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Fig. 58 Harmonic Analysis of the Sweep Effect on the Moment Response at =9 deg, 

k c = 0.075 (f =8 Hz, M c = 0.40), and a = 8 deg 
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